Collective flow of open and hidden charm in Au+Au collisions at y/s = 200 GeV 
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We study the collective flow of open charm mesons and charmonia in Au+Au collisions at ^/s — 
200 GeV within the hadron-string-dynamics (HSD) transport approach. The detailed studies show 
that the coupling of D, D mesons to the light hadrons leads to comparable directed and elliptic 
flow as for the light mesons. This also holds approximately for J/^ mesons since more than 50% 
of the final charmonia for central and mid-central collisions stem from D + D induced reactions in 
the transport calculations. The transverse momentum spectra of D, D mesons and J/^'s are only 
very moderately changed by the (pre-)hadronic interactions in HSD which can be traced back to 
the collective flow generated by elastic interactions with the light hadrons. 

PACS numbers: 25.75.-q, 13.60.Le, 14.40.Lb, 14.65.Dw 



I. INTRODUCTION 

The dynamics of nucleus-nucleus collisions at 
Relativistic-Heavy-Ion-Collider (RHIC) energies are of 
fundamental interest with respect to the properties 
of hadronic/partonic systems at high energy densities. 
Especially the formation of a quark-gluon plasma (QGP) 
and its transition to interacting hadronic matter has 
motivated a large community for more than two decades 
d S 01 ■ However, the complexity of the dynamics has 
not been unraveled and the evidence for the formation 
of a QGP and/or the properties of the phase transition 
is much debated Apart from the light and strange 
flavor (u,u,d,d,s,s) quark physics and their hadronic 
bound states in the vacuum {p,n,Tr,K,<f>,A etc.), the 
interest in hadrons with charm (c, c) has been rising 
continuously since the heavy charm quark provides an 
additional energy scale, which is large compared to 
Aqcd- Accordingly, the hadronic bound states (with 
a c or c quark) have a much larger mass than the 
ordinary hadrons and one has speculated that especially 
charmonia (cc bound states) might only be formed in 
the very early phase of the heavy-ion collision. 

In the past the charmonia J/^, Xo have been dis- 
cussed in context with the phase transition to the QGP 
since cc states might no longer be formed due to color 
screening |6j,|7j . However, more recent calculations within 
lattice QCD (LQCD) have shown that at least the J/ 1 ^ 
survives up to ~ 1.5 T c (T c sa 0.17 GeV) such that the 
lowest cc states remain bound up to energy densities of 
about 5 GeV/fm 3 It is presently not clear if 

also the D or D* mesons will survive at temperatures 
above T c but strong correlations between a light quark 
(antiquark) and a charm antiquark (quark) are likely to 
persist also above T c . 

Moreover, it has been pointed out (within statisti- 
cal models) that at top RHIC energies the charmo- 
nium formation from open charm and anticharm mesons 
might become essential 0, 0, 0, \^ an ^ even exceed 
the yield from primary nucleon-nucleon (TV TV) collisions 



|l5| . Such concepts should work out if the early hot 
and dense medium created in nucleus-nucleus collisions 
is very strongly interacting and an approximate chemi- 
cal equilibrium is achieved rapidly. As argued in Refs. 
16, 17, 181 an early equilibration might be due to multi- 
particle interactions either on the partonic |18| or the 
hadronic side |16lll7| |. 

A previous analysis within the HSD transport model 
|l9j | has demonstrated that the charmonium production 
from open charm and anticharm mesons becomes essen- 
tial in central Au+Au collisions at RHIC. This is in accor- 
dance with independent studies in Refs. |2jj,[2l|. On the 
other hand, these backward channels - relative to char- 
monium dissociation with comoving mesons - have been 
found to be practically negligible at SPS energies. Fur- 
thermore, the transport studies in Ref. |22| have shown 
that chemical equilibrium between the different charmo- 
nia J/^f, Xc is not obtained in full phase space on 
the basis of (pre-)hadronic interactions. As pointed out 
in the latter study this opens up a possibility to distin- 
guish experimentally a statistical freeze-out concept from 
a hadron/string dynamical picture. 

Apart from the total and relative abundancies of char- 
monia and open charm mesons also their collective prop- 
erties are of interest. Here especially the transverse mo- 
mentum (or mass) spectra are expected to provide valu- 
able insight to the dy namics in either the very early or 
late phase [H |H Hi |H H3 . We recall that the trans- 
verse mass spectra at midrapidity show experimentally 
a 'slope' parameter that increases (~ linearly) with the 
hadron mass [2^ if the latter involve only light quarks 
or only a single strange quark. Multi-strange baryons 
indeed show a 'slope' parameter below the linear scaling 
with mass known from the lighther hadrons |29j . The 
question thus arises if the open charm mesons will follow 
the trend of the light hadrons similar to kaons (involving 
only a single s or s quark)? Also: will the charmonia with 
their substantially higher mass also show the linear trend 
as expected from hydrodynamics [23 or do they decouple 
early as expected in hadron/string dynamical approaches 
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due to the rather low cross sections with light hadrons 
|3 (cf. also Ref. H)? 

In this work we extend our previous studies in Refs. 
[l9l I22I |27| with respect to the collective dynamics of 
D, D mesons and charmonia and concentrate in particu- 
lar on the transverse mass spectra, the in-plane flow v\ 
and the elliptic flow V2 of these particles. The question 
we aim at solving is: what is the amount of collectivity 
generated by (pre-)hadronic interactions in the course of 
the expansion of the system? Any sincere deviation to 
future measurements thus will indicate additional sources 
for pressure and/or strong interactions beyond the stan- 
dard hadron/string picture. 

Our study is organized as follows: After a very brief 
reminder of the conceptional organization of the HSD 
transport approach in Section II we will present the in- 
put for the transport calculations with respect to the 
open charm and charmonium dynamics (or provide the 
relevent references for details). In Section III we present 
our results for the transverse momentum spectra of D- 
mesons and J/^'s for different centralities in Au+Au 
collisions at y/s = 200 GeV. In order to see the effect 
of final state interactions more clearly we will, further- 
more, present ratios of momentum spectra from Au+Au 
reactions relative to scaled pp collisions. Section IV is de- 
voted to the in-plane flow v\ and elliptic flow V2 of open 
charm mesons and charmonia as well as to the freeze-out 
properties of these particles. Section V concludes this 
study with a summary. 



II. BASIC CONCEPTS OF THE HSD 
TRANSPORT APPROACH 

We employ the HSD transport model 0, H3 for our 
study of Au+Au collisions. This approach takes into ac- 
count the formation and multiple rescattering of formed 
hadrons as well as unformed 'leading' pre-hadrons and 
thus incorporates the dominant final state interactions. 
In the transport approach nucleons, A's, N*(1440), 
N*(1535), A, S and E* hyperons, S's, S*'s and O's as 
well as their antiparticles are included on the baryonic 
side whereas the 0~ and 1~ octet states are included in 
the mesonic sector. Inelastic hadron-hadron collisions 
with energies above Jls ~ 2.6 GeV are described by 
the FRITIOF model ^ (employing PYTHIA v5.5 with 
JETSET v7.3 for the production and fragmentation of 
jets [12), whereas low energy hadron-hadron collisions 
are modeled in line with experimental cross sections. We 
stress that no explicit parton cascading or gluon-gluon 
dynamics is involved in our transport calculations con- 
trary to e.g . the AMPT model [33 or explicit parton 
cascades |34| . 

A systematic analysis of HSD results and experimen- 
tal data for central nucleus-nucleus collisions from 2 
A -GeV to 21.3 A-TeV has shown that the spectra for 
the 'longitudinal' rapidity distribution of protons, pions, 
kaons, antikaons and hyperons are in reasonable agree- 



ment with available data [T^. I3H |3(|. However, there 
are problems with the dynamics in the direction trans- 
verse to the beam. Whereas the pion transverse momen- 
tum spectra are rather well described from lower AGS 
to top RHIC energies the transverse momentum slopes 
of kaons/antikaons are clearly underestimated above ~ 5 
A-GeV in central Au+Au collisions. In Ref. [3(j this fail- 
ure has been attributed to a lack of pressure generation 
in the very early phase of the heavy-ion collisions which 
also shows up in the underestimation of the elliptic flow 
of charged hadrons at RHIC energies ^3 • 

We note, that the inclusion of initial state Cronin ef- 
fects gives a substantial hardening of kaon spectra at 
RHIC energies, but the slope of the pion spectra at low 
tot is only slightly enhanced. In the present study we 
have also included the Cronin effect, however, found out 
that it has only a small impact on charm and charmonia 
spectra. 

Inspite of the deficiences pointed out above the overall 
reproduction of the experimental hadron spectra is suffi- 
ciently realistic such that we can proceed with open and 
hidden charm dynamics. 

A. Perturbative treatment of D, D mesons and 
charmonia 

The initial conditions for the production and also sub- 
sequent propagation of D,D mesons and charmonia are 
incorporated in the HSD approach by a superposition 
of pp collisions described via scaled PYTHIA [32| sim- 
ulations . For the production and propagation of 
open and hidden charm hadrons we employ a pertur- 
bative scheme as also used in Refs. 0, [22, H3- Each 
'perturbative' particle (hi) is produced in the transport 
calculation with a weight Wj given by the ratio of the 
actual production cross section divided by the inelastic 
nucleon-nucleon cross section, e.g. 

a NN ^ hi+x (s) 

VVl _inelas.f«,\ ' ^ 

a NN \ b ) 

We then follow the motion of the 'perturbative' parti- 
cles within the full background of strings/hadrons by 
propagating them as free hadrons, i.e. neglecting in- 
medium potentials, but compute their collisional history 
with baryons and mesons or quarks and diquarks. The 
actual parametrizations of the total cross sections for 
D, D* , D s , D* mesons as well as their antiparticles are 
given m Refs. 0H3 and [13 together with the cross 
sections for charmonium production in nucleon-nucleon 
and pion-nucleon collisions. The only modification intro- 
duced here - relative to Refs. 0, |53 - is a change of 
the parametrisation for the differential production cross 
section in transverse momentum. Our novel parametrisa- 
tion follows the first results from STAR [37j and is given 
by a power law in transverse momentum px 

W^T)~(l + f§)" 8 ' 3 (2) 



For simplicity - and lack of further information - we as- 
sume that this parametrization holds for all D, D mesons 
as well as charmonia. 



B. Cross sections involving charm hadrons 

In order to study the effect of rescattering we tenta- 
tively adopt the following dissociation cross sections of 
charmonia with baryons independent on the energy (in 
line with Refs. G3|H|27|): 

a cE B = 6 mb; Vj/^b = 4 mb; (3) 

<7 Xc b — 5 mb; aq,>B = 5 mb. 

In (JjJJ the cross section <j C cB stands for a (color dipole) 
pre-resonance (cc) - baryon cross section, since the cc 
pair produced initially cannot be identified with a par- 
ticular hadron due to the uncertainty relation in energy 
and time. For the lifetime of the pre-resonance cc pair 
(in it's rest frame) a value of r cg = 0.3 fm/c is assumed 
following Ref. j3g. This value corresponds to the mass 
difference of the and J . 

For D, D* , D, D* - meson (tt, r\, p. w) scattering we ad- 
dress to the calculations from Ref. [2]], [3t| which predict 
elastic cross sections in the range of 10-20 mb depending 
on the size of the formfactor employed. As a guideline 
we use a constant cross section of 10 mb for elastic scat- 
tering with mesons and also baryons, although the latter 
might be even higher for very low relative momenta. 

The cross sections for charmonium production by open 
charm mesons or the inverse comover dissociation cross 
sections are not well known and the significance of these 
channels is discussed controversely in the present litera- 
ture (cf. Ref. and Refs. therein). We here follow 
the simple 2-body transition model introduced in Refs. 
|l9ll22j | with a single free parameter, i.e. a matrix element 
squared Mq, that allows to implement the backward re- 
actions uniquely by employing detailed balance for each 
individual channel. The free matrix element has been 
fixed in Ref. (2^ in comparison to Jf^l and \P' suppres- 
sion data at SPS energies. For further details we have 
to refer the reader to Refs. 0, in order to avoid 
unnecessary repetitions. 



III. TRANSVERSE MOMENTUM SPECTRA 

Since the results for charmonium suppression at SPS 
and RHIC energies - within the present approach - have 
been presented in Refs. |l9ll22| we directly continue with 
the transverse dynamics of open charm mesons and char- 
monia. In Fig. ^we present the transverse px spectra of 
the final D + 5-mesons (solid lines with full dots, color: 
blue) and J /if? (solid lines with open squares, color: red) 
from Au + Au collisions at y/s — 200 GeV calculated 




p T [GeV/c] p T [GeV/c] 



FIG. 1: (Color online) The transverse pt spectra of the final 
D + _D-mesons (solid lines with full dots, color: blue) and 
J '/vp (solid lines with open squares, color: red) from Au + 
Au collisions at ^/s = 200 GeV calculated for b = 1, 5, 7 and 
12 fm at midrapidity. The dashed lines correspond to the pr 
spectra of D-mesons at the production point. 



for impact parameter b — 1,5,7 and 12 fm at midra- 
pidity. The dashed lines correspond to the pt spectra 
of Z?-mesons at the production point for reference, i.e. 
without any final state interactions. For very peripheral 
Au+Au collisions we find no effects from final state in- 
teractions whereas with increasing centrality there is an 
enhancement of the spectra at low and moderate pt from 
Au+Au collisions as well as a modest suppression at high 
Pt- 

We note that the power low parametrization gives 
a substantial enhancement of the high momentum tail of 
the open charm and charmonia px spectra compared to 
the previous parametrization used in HSD |l9l |22| . This 
is an improvement in the present approach as the new 
distribution provides a higher number of charm pairs at 
the beginning of the collisions. This will also affect the 
J /if? yield. Further experimental data on high px spectra 
in elementary as well as heavy-ion collisions will fix the 
proper shape. 

In order to quantify the effect of final state interactions 
we show in Fig. [21 the ratio of the final to the initial (i.e. 
at the production point) transverse px spectra of D + D- 
mesons (solid lines with full dots, color: blue) and J /if? 
(solid lines with open squares, color: red) from the same 
reaction for impact parameter b — 1,5,7 and 12 fm at 
midrapidity. 

The actual ratios in Fig. [3 show an enhancement of 
D, D mesons at low momenta with a maximum at px ~ 
1 GeV/c and a relative suppression for px > 2 GeV/c. 
These effects increase with the centrality of the Au+Au 
collision. For the J/\& a maximum in the ratio shows up 
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FIG. 2: (Color online) The ratio of the final to the initial 
(i.e. at the production point) transverse pt spectra of D + D- 
mesons (solid lines with full dots, color: blue) and J/^ (solid 
lines with open squares, color: red) from Au + Au collisions 
at y/s = 200 GeV calculated for 6 = 1,5,7 and 12 fm at 
midrapidity. 



at ~ 2 GeV/c while it drops below 1 forpr > 3-4 GeV/c. 

In accordance with the ratios in Fig. 2 the averaged 
value of the transverse momentum px, which is displayed 
in Fig. |3 at midrapidity for D + D-mesons (solid line 
with full dots, color: blue) and J/^f (solid line with 
open squares, color: red) from Au + Au collisions at 
y/s = 200 GeV, slightly increases (decreases) for J/^J> 
(D, D mesons) with decreasing impact parameter b. The 
-D-meson spectrum is relatively 'hard' due to the large in- 
trinsic charm quark mass. Through scatterings, the high 
Pt D-mesons loose their energy and develop the collec- 
tive shoulder at low px- This explains the dropping of 
D-meson (pt) in more central Au + Au collisions shown 
in Fig. |21 Such hadronic energy loss for .D-mesons will 
complicate the analysis of the gluon radiative energy loss 
before hadronization |4jj. These main trends should be 
readily tested experimentally in the near future. 

The question comes up how to interpret the findings 
in Figs, n to [31 We note that the maxima in the D + D 
and J/ty ratios disappear when switching off the rescat- 
tering with mesons in the transport approach. Thus a 
collective acceleration of the D + D mesons occurs domi- 
nantly via elastic scattering with mesons. This argument 
does not hold at first sight for the J/^F's since their cross 
section with mesons is much lower and no strong col- 
lective acceleration by rescattering should be expected 
|27|. The answer comes about as follows: As known from 
Refs. [l9l |22| the major fraction of final J/^P's stems 
from D + D — ► J/^(Xc) + meson channels and not from 
the primary production by NN collisions. This finding is 
quantified in Fig. 0] where the channel decomposition for 
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FIG. 3: (Color online) The average value of the transverse 
momentum pr at midrapidity for D, Z)-mesons (solid line with 
full dots, color: blue) and J/ty (solid line with open squares, 
color: red) from Au + Au collisions at yfs = 200 GeV versus 
impact parameter 6. 



the final J/^'s is shown as a function of the impact pa- 
rameter b in Au+Au collisions at y/s — 200 GeV. Except 
for peripheral reactions (b > 9 fm) the D + D channel 
(solid line with full dots, color: red) dominates while the 
production in initial baryon-baryon (BB) scattering (full 
line with open squares, colonblue) superceeds all other 
reaction channels for the most peripheral reactions. The 
contribution from meson-baryon (mB) channels (solid 
line with full triangles, color: green) is seen to be low 
at all centralities. 
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FIG. 4: (Color online) The channel decomposition for the fi- 
nal J '/SP in Au + Au collisions at yfs — 200 GeV within the 
transport approach as a function of the impact parameter b. 
The fraction of the final J/ty mesons produced in BB inter- 
actions is shown by the solid line with open squares (color: 
blue), the contribution from meson-baryon (mB) interactions 
is given in terms of the solid line with full triangles (color: 
green) and the DD annihilation fraction is displayed by the 
solid line with full circles (color: red). 

We recall that the sensitivity of J/\V formation versus 
D + D annihilation to the formation amplitude has been 
studied in detail in Refs. [l9l I22]]. The amplitudes for 
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D + D < - > J/*(xc, + mesons have been fitted to 
the experimental data on J/^S suppression and iff' / Jty 
ratios at SPS energies and used also for the RHIC ener- 
gies. This is legitimate because the latter reactions are 
typical 'comover' reactions with low/moderate relative 
momenta. The average relative momenta for comover re- 
actions do not change very much from top SPS to RHIC 
energies. At RHIC energies, however, the density of pro- 
duced D-mesons is much higher ( 17 D + D pairs in cen- 
tral collisions), and thus the probability to annihilate to 
charmonia (~ p 2 ) is dramatically higher (cf. Fig. 6 in 
Ref. |22j). In fact, the system is close to chemical equi- 
librium in central reactions such that a sensitivity to the 
matrix elements is lost. 

According to Fig. 21 the bulk part of the J/^f's thus 
stems from D + D induced channels and their produc- 
tion happens somewhat delayed in time when the D, D 
mesons have already picked up collective transverse mo- 
mentum. In the D + D — ► J/^{Xc) + meson reaction 
then the formed J/$ carries the collective momentum of 
both D and D such that the charmonium appears to be 
accelerated even more than the £>-mesons. On the other 
hand the J/'&'s with a high transverse momentum (in 
the power law tail) partly get dissociated with baryons 
of the target /projectile or with mesons in the expansion 
phase of the system. The probability to recreate a high 
px charmonium by the D + D channel is lower than in 
the initial TV TV channel since the invariant energy y/s in 
the D + D reaction is lower by more than an order of 
magnitude and exponential (thermal) J spectra be- 
come populated. This interpretation will, furthermore, 
be supported by the studies in the next Section. 



IV. DIRECTED AND ELLIPTIC FLOW 

Apart from the transverse flow that shows up as a 
'shoulder' in the px spectra at low momentum or a max- 
imum in the ratio relative to pp spectra (cf. previous 
Section) the directed flow 



Vi(v,Pt) 
and the elliptic flow 

V2(y,Pr) = 



Px\ 




(4) 



(5) 



V,Pt 



provide additional information on the collective currents 
of hadrons in the complex reaction |41| . 

Fig. [SJshows the HSD predictions for the directed flow 
v\ (upper part) and elliptic flow v-i (lower part) of D + D- 
mesons (solid lines with full dots) and J/^f (solid lines 
with open squares) from Au + Au collisions at y/s = 200 
GeV for b — 7 fm versus px for < y < 1 (left part) and 
rapidity y (right part) integrated over px- 

The directed flow v\ is negative (within statistics up to 
3.5 GeV/c) for both D + D mesons and J/^'s in the 
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FIG. 5: (Color online) The directed flow v\ (upper part) and 
elliptic flow V2 (lower part) of D + D-mesons (solid lines with 
full dots, color: blue) and J/SP (solid lines with open squares, 
color: red) from Au + Au collisions at ^/s = 200 GeV for 
b = 7 fm versus pr for < y < 1 (left part) and rapidity y 
(right part) integrated over pr- 



rapidity interval < y < 1. There is a tendency that v± 
for J/'f's is larger in magnitude for 0.5 GeV/c < px < 2 
GeV/c which supports the interpretation given in Section 
III that the final J/^'s dominantly stem from D + D 
channels where both open charm mesons already have 
picked up some collective flow. The pronounced 'antiflow' 
of open charm mesons and J/vP's becomes visible in the 
r.h.s. of Fig. [5] (upper part) where the slope dv\/dy at 
midrapidity for J/'J's is also slightly larger in magnitude 
than that for D + D mesons. 

The question now emerges how the 'bulk' of the light 
hadrons flows in this reaction. The answer is given in 
Fig. [S] where the HSD result for the directed flow v\ for 
charged hadrons from semi-central Au + Au collisions at 
y/s = 200 GeV is plotted versus pseudorapidity n and 
compared to the data from the STAR Collaboration j3?J 
(solid dots) and the PHOBOS Collaboration jH (solid 
triangles). Apparently the charged hadrons show a com- 
parable flow Vi (rj) as the open charm mesons and char- 
monia; consequently, the latter flow with the bulk of the 
lighter hadrons made up from u, d, s and u, d, s quarks. 
We recall that the HSD calculations for the directed flow 
of charged hadrons match reasonably well with the data 
from Refs. |H for \n\ < 3 (cf. Ref. @ for related 
UrQMD calculations). 

Coming back to the results for the elliptic flow V2 in the 
lower part of Fig. [3]we notice that both J/\l/'s and D, D 
mesons show in-plane flow since V2(y,PT) > 0. The ellip- 
tic flow of D,D mesons is larger for 1 GeV/c < px < 3 
GeV/c than that for J/^'s which indicates that the D, D 
mesons are accelerated earlier than the J/^'s. The el- 
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FIG. 6: (Color online) The directed flow Wj for charged 
hadrons from semi-central Au + Au collisions at ^J~s — 200 
GeV versus pseudorapidity u in comparison to the data from 
STAR :35 (solid dots) and PHOBOS (solid triangles). 



liptic flow V2 shows a maximum around midrapity as the 
charged hadrons [I^ and is also slightly larger for D, D 
mesons than for J/\I/'s. However, for impact parameter 
b = 7 fm the elliptic flow of open charm mesons is < 3% 
whereas the elliptic flow of charged hadrons reaches up 
to ~ 5% at midrapidity. We recall that the HSD calcula- 
tions underpredict the vi of charged hadrons at midrapid- 
ity by about 30-35% 0. Consequently one should ex- 
pect also a larger elliptic flow for the open charm mesons 
and charmonia in experiment. We mention that in the 
quark coalescence model of Ref. [25| a significantly larger 
elliptic flow is obtained for both D,D mesons and J/^>. 

We now turn to the multiplicity of D, D mesons and 
J/4" versus centrality in Au+Au collisions at ^fs = 200 
GeV. The left part of Fig. [7| shows the multiplicity of 
D + D pairs (solid lines with full dots, color: red) and 
J/vE' (solid lines with open squares, color: blue) from 
Au + Au collisions at y/s — 200 GeV integrated over 
rapidity versus impact parameter b, whereas the right 
part displays the ratio of the multiplicities of J/>F over 
D + D pairs versus b. Note that the J/'B multiplicity on 
the l.h.s. has been multiplied by a factor of 100. Whereas 
we expect ~ 17 D + D pairs in the top central collisons 
the J/'F multiplicity is only ~ 5xl0~ 2 , i.e. lower by a 
factor of ~ 300. 

The r.h.s. of Fig. shows the ratio of the multiplicity 
of final J/vP's over D + D pairs versus impact parameter 
b. Since the multiplicity of D + D pairs scales directly 
with the number of binary (initial) AW collisions in the 
HSD transport approach |27| this ratio gives an infor- 
mation about J /\& suppression relative to binary scaled 
pp collisions. We find - in accordance with the previous 
studies in Rcfs. |19l l22l | - that the J/'F suppression in 
very central collisions relative to very peripheral reactions 
is about a factor of ~ 2.5. Note that the explicit shape 
of this ratio versus b is the result of a complex coupled- 
channel problem that cannot be anticipated by simple 
scaling arguments. We recall that the dominant charmo- 
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FIG. 7: (Color online) Left part: the multiplicity of D + D 
pairs (solid lines with full dots, color: red) and J/$ (xlOO, 
solid lines with open squares, color: blue) from Au + Au 
collisions at i/s = 200 GeV integrated over rapidity versus 
impact parameter b. Right part: the ratio of the multiplicity 
of J/$ over D + D pairs versus b. 



nium absorption channel in our transport calculations 
is the very early color-dipole dissociation with nucleons 
of the target/projectile (or their baryon-like remnants). 
These early dissociation reactions involve a high invari- 
ant collision energy such that the final fragments are dis- 
tributed in a wide rapidity range which makes any recom- 
bination back to the original color-dipole - nucleon chan- 
nel very unlikely. The formation of charmonium states 
(as well as D, D mesons ) occurs later with characteristic 
formation times of 0.3 and 0.5 fm/c in their rest frame, 
respectively. Moreover, any dissociation with mesons is 
delayed until the latter have formed (cf. ^H). Now the 
dissociation of charmonia with comoving hadrons is com- 
pensated to a large extent by the inverse D + D — > char- 
monium + meson channels at top RHIC energies. The 
actual time integrated reaction rates for the forward and 
backward channels are shown in Fig. 5 of (2^ versus 
centrality. Since both reaction channels turn out to be 
comparable in the time integrated rate for midcentral and 
central collisions, the net result is not evident a priori. 
We recall that the J/'F's gain slightly from the D + D 
channels whereas the \c an d ^' states loose from the 
dissociation channels with mesons (cf. Fig. 6 in H2). 
Now the feeddown from \c and '5' to J/'F's has to be in- 
cluded for the final J/ty multiplicity which for Au + Au 
at top RHIC energies shows a small net suppression from 
comover reactions on top of the early suppression by in- 
teractions with projectile/target nucleons. This rather 
complicated coupled-channel scenario finally leads to the 
ratio of J/'F over D + D pairs in Fig. [7| 

We now turn to the freeze-out properties of D + D 
mesons and J/'F's. The upper part of Fig. [S] presents 
the distribution in time of the last interaction for the final 
J/4" mesons (xlOO, solid line, color: blue) and D + D 
mesons (dashed line, color: red) from Au + Au collisions 
for 6 = 2 fm at V« = 200 GeV at midrapidity (\y\ < 1). 
Here the production time is also recorded if the particles 
did not suffer from a further interaction (initial peaks up 
to 2 fm/c). The time t = is defined by the contact time 
of the heavy ions. Whereas the J/^S mesons freeze-out 



7 



0.5 



1 04 - 



5 0.3 



a 

? 0.2 



0.1 _ 



z 

T3 



0.0 



800 



600 



* 400 



200 







~i 1 1 1 1 1 1 1 1 1 1 1 r 

Au+Au, s" 2 =200 GeV, b=2 fm, lykl 




■finalj/* * 100 
• final D+Dbar 



_L. 



J_ 







25 



30 



10 15 20 25 30 35 
time [fm/c] 







1 1 1 1 1 
Au+Au, 


i 1 i 1 i 1 i 1 

s" 2 =200 GeV, b=2 fm, lykl - 
Collision rate of 






M 

^ lAtotal 


formed + leading particles 


















mm > 

A 














V 


»«">.-*, 


- j - - : i i — j 



35 



5 10 15 20 
time [fm/c] 

FIG. 8: (Color online) Upper part: the distribution in time 
of the last interaction for the final J/tp mesons (xlOO, solid 
line, color: blue) and D + D mesons (dashed line, color: red) 
from Au + Au collisions for b = 2 fm at y/s — 200 GeV (\y\ < 
1). Here the production time is also recorded if the particles 
did not suffer from a further interaction (initial peaks up to 
2 fm/c). Lower part: the time evolution of the collision rate 
of formed and leading particles from baryon-baryon (BB), 
meson-baryon (mB) and meson-meson (mm) collisions from 
Au + Au collisions for b = 2 fm at — 200 GeV integrated 
over rapidity. The time t = is denned by the contact time 
of the heavy ions. 



a quark from a struck nucleon - forming an endpoint of 
the excited string - interacts without time delay (cf. Ref. 
|44l | for a detailed discussion of the pre-hadron concept 
in HSD). Furthermore, the early interactions of formed 
mesons (for 1 fm/c < t < 4 fm/c) occur dominantly in 
the outer transverse region of the initial fireball where 
the energy density is below 1 GeV/fm 3 . 

As seen from Fig. 8 the mm collision rate of formed 
mesons also starts with a short delay in the order of the 
formation time Tf but superceeds the mB collision rate 
by more than a factor of 4. The maximum in the mm 
collision rate occurs for t w 3 — 4 fm/c when the energy 
density in a sizeable volume drops below 1 GeV/fm 3 and 
the system is still very dense. We point out, further- 
more, that the 'hadronic burning' by meson-baryon and 
meson-meson collisions persists up to rather long times 
since bunches of close-by hadrons may interact almost 
continuously if their relative momenta are low. 

When comparing the upper part of Fig. |H1 with its 
lower part one notices that the maximum in the mm col- 
lision rate occurs much earlier than the average freeze- 
out of D + D mesons (at 9-12 fm/c). This further sup- 
ports our conjecture in the end of Section III that the 
D+D mesons suffer a couple of interactions with the light 
mesons and pick up collective flow before they freeze-out 
rather late on a time scale that is only slightly shorter 
(but comparable) to the freeze-out time for light hadrons 
if one discards the very late decay of resonances (e.g. 
u) -> 3tt's, A -> Nn, K* ->■ Kir etc.). 



V. SUMMARY 

In this work we have extended our previous inves- 
tigations within the HSD transport approach in Refs. 
fl9l I22I |27| with respect to the collective dynamics of 
D, D mesons and charmonia. Our detailed studies of 
Au+Au collisions at y/s = 200 GeV have shown that: 



early due to their small cross section with light hadrons 
the D + D mesons interact for a long time and freeze-out 
on average at 10 — 12 fm/c. 

The lower part of Fig. [S] shows the collision rate of 
formed and leading particles from baryon-baryon (BB), 
meson-baryon (mB) and meson-meson (mm) collisions 
from Au + Au collisions at b = 2 fm and \fs = 200 
GeV integrated over rapidity. Here leading diquarks are 
counted as baryons (B) whereas leading quarks and an- 
tiquarks are registered as mesons (m). The rate of BB 
collisions shows a strong peak during the passage time 
of the heavy ions (< 0.2 fm/c) and is almost negligible 
later on. The mB collision rate (apart from the initial 
peak for t < 0.2 fm/c) shows a sizeable contribution only 
after a short time delay of ~ 1 fm/c which is the scale 
of the hadron formation time for low momenta in the 
cms. We point out that the very early mB collisions (t < 
0.4 fm/c) correspond to pre-hadron interactions, where 



• Pre-hadronic and hadronic interactions generate a 
transverse collective flow of D, D mesons and char- 
monia that shows up as a 'shoulder' in the low mo- 
mentum spectra or as a maximum in the ratio rel- 
ative to scaled pp collisions. 

• The high p^ power law tail of the spectra is only 
moderately suppressed. 

• The directed flow of D, D mesons and charmonia is 
comparable to that of the light charged hadrons. 

• The elliptic flow of D,D mesons and charmonia is 
smaller than that of the light hadrons. 

• The D, D mesons freeze out on average at ~ 9 — 12 
fm/c which is shorter than the freeze-out time for 
light hadrons (when neglecting explicit resonance 
decays). 
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Whereas the collective acceleration of the D + D 
mesons via elastic scattering with mesons in the ex- 
pansion phase of the 'fireball' can be well understood 
on the basis of the rather large cross section of D, D 
mesons with nucleons and mesons |2ll l39j - in view of 
the light quark/ antiquark content of these states - the 
collective dynamics of the J/'f's cannot be explained 
in this way since the J/^> cross section with mesons is 
much lower. However, the dominant fraction ( > 50%) 
of final J/^'s in central and mid-central reactions stems 
from D + D — ► J/^{Xc) + meson channels and to a 
lower extent from the primary production by NN col- 
lisions (cf. Fig. 4). Thus the production of the fi- 
nal charmonia happens delayed in time when the D, D 
mesons have already picked up collective flow from in- 
teractions with light (and fast) hadrons. Furthermore, 
in the D + D — > J/^>(x c ) + meson reaction the formed 
J/vP carries the collective momentum of both D and D 
such that the charmonium appears to be accelerated even 
more than the D- mesons (cf. Figs. 2 and 5). On the 
other hand the elliptic flow of D,D mesons is slightly 
larger than that of J/^f's which indicates on average an 
earlier production of the open charm mesons. 



The differential spectra, ratios as well as the dif- 
ferential flow analysis for Ui(j/,J>t) an d V2(y, pr) for 
open charm hadrons and charmonia is expected to be 
proved/disproved in the near future by the RHIC exper- 
iments. Our analysis and interpretation of the results 
has also paved the way to extract relevant time scales 
and possibly interaction rates or freeze-out times. Sen- 
sible deviations from our predictions will points towards 
a dynamical origin that is not included in our present 
calculations and should be addressed to explicit partonic 
interactions in a possibly colored medium. 
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